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Abstract 

Combining the WKB expansion at large distances and Perturbation Theory at small distances 
it is constructed a compact uniform approximation for eigenfunctions. For lowest states Isag 
and 2pau this approximation provides the relative accuracy < 10^^ (5 s.d.) for any real x in 
eigenfunctions and for total energy E{R) it gives 10-11 s.d. for internuclear distances R G [0,50]. 
Corrections to proposed approximations are evaluated. Separation constants and the oscillator 
strength for the transition IsGg — )• 2pau are calculated and compared with existing data. 
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INTRODUCTION 



The molecular ion is the simplest molecular system which exists in Nature. It was 
the first studied molecular system since the inception of the quantum mechanics which later 
appears in all QM textbooks (see e.g. ll). Needless to say that this system plays very 
important role in different physical sciences, in particular, in laser and plasma physics. 

From technical point of view this is the unique molecular system which admits complete 
separation of variables (in elliptic coordinates). Definitely, this problem is non-solvable. 
Thus, the problem can be solved in approximate way only. We introduce a natural definition 
of solvability of non-solvable spectral problem: for any eigenfunction \E' we can indicate 
constructively an uniform approximation "^app such that 



< ^ , (1) 



in the coordinate space. It implies that any observable, any matrix element can be found with 
accuracy not less than 6. A simple idea we are going to employ is to combine WKB-expansion 
at large distances with perturbation theory at small distances near extremum the potential 
in one interpolation. Recently, this idea was realized for quartic anharmonic oscillator ^ 
and double- well potential jsj. In both cases for the lowest states it was constructed two- 
three parametric uniform approximations of the eigenfunction leading to 10 s.d. in energies 
and with 6 ~ 10~^ for any value of the coupling constant and size of the barrier. The goal 
of this Letter is to present such an approximation with 6 ~ 10~^ for two lowest (and the 
most important) states Isag and 2pau of the molecular ion. It is worth mentioning 
that a study of the wavefunctions of the molecular ion in a form of expansion in some 
)asis was initiated by Hylleraas j4| and was successfully realized in the remarkable paper 
(see also 0, B])- Attempts to find bases leading to fast convergence are still continuing. 
At present, the pure exponential basis seems the most fast convergent (see e.g. 1^8,] and 
references therein). It is worth noting that following the analysis of classical mechanics of 
the system and its subsequent semiclassical quantization some uniform approximations 

n 

of wavefunctions of low lying electronic states were constructed [9] . Local accuracies of these 
approximations are unclear albeit eigenparameters are found with a few significant figures. 
The Schrodinger equation, which describes the electron in the field of two centers of the 
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charge Z = 1 at the distance R, is of the form 

A - — - — ) ^ = ^'^ , ^ G L^(R^) , (2) 
n r2 J 

where E' = [E — -|) and the total energy E are in Rydbergs, ri 2 are the distances from 
electron to first (second) center, respectively. Following [it] let us introduce the dimensionless 
elliptic coordinates [l^: 

i = ^ , V = ^ , 1 < ? < 00 , -l<r]<l (3) 

and azimuthal angle (j). The Jacobian is oc {^"^ — rf). The equation ([2]) admits separation 
of variables in ([3]). Since the pro ject ion of the angular momentum to the molecular axis 
commutes with the Hamiltonian [iSj the eigenstate has a definite magnetic quantum number 
A. The Hamiltonian is permutationally-symmetric ri ^ r2, or, equivalently, t] — )■ —77, hence, 
any eigenfunction is of a definite parity (±). As a result, it can be represented in a form 

* = X{i){e - lt'^Y{r^){l - rf)^l\^^^t> , A = 0, 1,2, . . . (4) 

where Y{r]) is of definite parity. After substitution of (j4]) into (j2]) we arrive at the equations 
for X(0 and Y{r]), 

- l)d^X] + 2Kid^X + ^2R^ + A]X = , Xe L\C e [1, 00)) , (5) 

d^[{7]^ - l)d.,Y] + 2A7]dr,Y + [-pY + A]Y = 0,Ye L^t] G [-1, 1]) , (6) 
respectively, where following js] we denote, 

P = ^ , (7) 

and A is a separation constant. Equations (jS]), dHD define a bispectral problem with E, A as 
spectral parameters. Both spectral parameters E, A depend on R. Square-integrability of 
the function \E' (jlj) implies a non-singular behavior of X at — )■ 1 as well as non-singular 
behavior of y at 77 —t- ±1. The latter condition implies a certain behavior of the solution 
Y at large arguments 77. A non-singular solution 1^(77) can be unambiguously continued in 
77 beyond the interval [—1,1], it has to be growing (non-decaying) at |r7| 00. It is in 
agreement with large-?] behavior of the Hund-MuUiken function (it mimics the incoherent 
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interaction of electron with charged centers) for both Isag (parity +) and 2pau (parity -) 
states 

XT/Ci) _ p-2a2ri i -2a2r2 _ n -a2R( COsh(o;2-R^) 
^ HM — ^ =■= *^ — 

sinh(a2-R'7) 

which describes large R behavior, similarly, for the Guillemin-Zener function (it mimics the 
coherent interaction of electron with charged centers) 



(8) 



(±) 

GZ 



g-2o3ri-2o4r2 _|_ g-2o3r2-2a4ri _ 2e~('^3+°4)fiC 



cosh((a3 — a^jRr]) 
sinh((Q;3 — ais)Rr]) 



(9) 



which has to correspond to small R behavior. 
Asymptotics. If we put X = e"*^, then at — )■ oo. 



^ = p^-(^-A-i) ioge + 



^ + (f-A-l)(f + A) 



P 



P 



1 



(10) 



which is nothing but WKB-expansion, and at ,^ — )■ 0, 



2^ 3^ 



(p2 + A^- A{2A + 3)) 



12 



(11) 



Similarly to X we put Y = e ^, then at 77 — )■ 00, 



,'A-A{A + 1) \ 1 
g = -pv + (A + 1) log ?7 - ( p\ ^ + • • • 



when at 7] — )■ 0, 



A 2 , {p' + A'-A{2A + 3)) , 

V H r) + . . . . 

2 ' 12 ' 



(12) 



(13) 



The important property of the expansions ( jlOl) and (1121) is that the coefficients in front of 
the growing terms at large distances (linear and logarithmic) are found explicitly, since they 
do not depend on the separation constant A. 

Approximation. Making interpolation between WKB-expansion (fTOj) and the perturbation 
theory flTTl) for X, f[T2]) and f[T^ for F, correspondingly, and taking into account that the 
Z2-symmetry of \E': rj ^ —i] is realized through use of cosh(sinh)-function (cf. ([H]) and (^) 
we arrive at the following expression 



(±) 

n,m,A 



ie - i)^/^P„.(0 {i-vY^'Qmin 

5 — e T+« 

(7 + 0^' 



l+A- 



;i + b2V^ + hrj^ 



1+A 



cosh / ai + pa2V + phri^ 
V 

sinh V 



1 + b2V^ + hr]^ 



(14) 
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for the eigenfunction of the state with the quantum numbers (n,m,A,±). Here 0,7 and 
oi,2, ^2,3 are parameters (see below), Pn{0 Qm{,if) are some polynomials of degrees n and 
m with real coefficients with n and m real roots in the intervals [1, 00) and [0, 1], respectively. 
These polynomials should be chosen in such a way to ensure their orthogonality. 
Results. As an illustration we consider two lowest states - one of positive and one of negative 
parity, Iscr^ (0,0,0,+) and 2pau (0,0,0,—), respectively. Corresponding approximations 
have the form 



^0,0 = « 



(7 + 0'-f (i + M' + M')^/' 



cosh f ai+ pa2rf + ph^r]'^ 



sinh V 1 + + M 
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^X,{OYt\v) , (15) 

(cf. ffT^ ) and each of them depends on six parameters 0,7 and ai^2,&2,3- The easiest way 
to find these parameters is to make a variational calculation taking (fT5|) as a trial function 
for R fixed and with p as an extra variational parameter. Immediate striking result of the 
variational study is that for all R G [1,50] the optimal value of the parameter p coincides 
with the exact value of p (see (jTj)) with extremely high accuracy for both Isag and 2pau 
states. It implies a very high quality of the trial function - the variational optimization wants 
to reproduce with very high accuracy a domain where the eigenfunction is exponentially 
small, hence, the domain which gives a very small contribution to the energy functional. 
In Tables I, II the results for the total energy (as well as for sensitive p) vs R of Iscr^ and 
2pau states are shown as well as their comparison with ones obtained by Montgomery in 



and also with the results we 



highly-accurate realization of the approach by Bates et al 
obtained in the Lagrange mesh method based on Vincke-Baye approach [1^ (details will be 
given elsewhere). For all studied values of R for both Iscr^ and 2pau states our variational 
energy turns out to be in agreement on the level of 10 s.d. with these two alternative 
calculations. Variational parameters are smooth slow-changing functions of i?, see Tables 
III-IV. All calculations were implemented in double precision arithmetics and checked in 
quadruple precision one. It is worth noting that the number of optimization parameters can 
be reduced putting 02 = 62 = - the accuracy in energy drops from 10-11 to 5-6 significant 
digits. 

Hence, our relatively-simple, few parametric functions f|T5|l taken as trial functions in 
a variational study provide very high accuracy in energy. The natural question to ask is 



how close these functions are to the exact ones in configuration space. In order to study 
this question we develop a perturbation theory in the Schroedinger equation ([2]) taking 
a trial function (IT^ as zero approximation. The easiest way to realize it is to consider 
non-linearization procedure 13|: if the potential is of the form V = Vq + XVi, then it is 
looked for energy and the eigenfunction in the form of power-like series in the parameter A, 
E = ^ A^i^n and \1/ = \l/o exp(— ^ A"'^?^), respectively. Due to specifics of (1) because of 
the separation of variables the procedure can be developed for both functions X and Y (see 
(jl])) separately as well as for the separation parameter A, while keeping the energy E fixed. 
It can be done for the system of equations (jSD, dn])- For simplicity we consider nodeless in ^ 
and 77 states, (0,0,A, ±). As a first step let us transform (jS]) into the Riccati form by 
introducing X = e~'^ and Y = e~^, respectively, 

{e - l){x' - x^) + 2(A + 1)C X = A- V{0 , x = ^'^ (16) 

where the "potential" V{^) = a^i^ - 2i?^, and 

{rf - l){y' - y') + 2{A + l)r] y = A - W{v) , y = g', (17) 

where the "potential" Wlr]) = a^Tf. 

Let us choose some xo(0 = V'olO' then substitute it to the l.h.s. of flT6|) and call the result 
as unperturbed "potential" Vo(0 putting without loss of generality Aq = 0. The difference 
between the original V{C,) and generated Vo(0 is the perturbation, Vi{^) = V{^) — Vo{^). For 
a sake of convenience we can insert a parameter A in front of Vi and develop the perturbation 
theory in powers of it, 

x = 5^A"x„ , A = 5^A"A„,, . (18) 
The equation for nth correction has a form, 

ie - 1)< - 2[{e - l)xo - (A + l)C]xn = - (19) 
where Vi = Vi and f„ = ((^^ — 1) Y17=i ^iXn-i for n > 1. It can be immediately solved, 

ie - i)^+^e-2^°x„ = ^Vn,. - v^){e - i)^e-2^°rfe , (20) 



and 



In a similar manner by choosing ?/o(^) = Qoiv)^ building the unperturbed "potential" Wo(ry) 
and putting = as zero approximation one can develop perturbation theory in the 
equation (ITTl) 

?/ = 5^A"2/„ , A = 5^A"yl„,, . (22) 
The equation for nth correction has a form similar to f|T9|) . 



ir]^-l)y'^-2[{7]^-l)yo-iA + l)v]yn = An,y - (23) 
where wi = Wi = W — Wq and Wn = (?7^ — 1) J2i=i UiVn-i for n > 1. Its solution is given by 
(^2 _ i)A+ig-2,o^^ = l\An,y-Wn){v'-lfe-'^'dr] , (24) 

(cf.d^nD) and 

-^n y ~~ J (^^) 

(cf. fl2T]) ). In order to realize this perturbation theory a condition of consistency should be 
imposed 

An,x — ■^n,y ■ (26) 

This condition allows us to find the parameter p and, hence, the energy E' and E (see ([7j)). 

Sufficient condition for such a perturbation theory to be convergent is to require a per- 
turbation "potential" to be bounded, 

|^i(OI <Q , |W^i(^)| <C, , (27) 

where C^, are constants. Obviously, that the rate of convergence gets faster with smaller 
values of C^.Cr^. It is evident that the perturbations Vi{^) and Wiirf) get bounded if v'o(^) 
and Qo{ri) are smooth functions vanishing at the origin but reproduce exactly the growing 
terms at \^\, \r]\ tending to infinity in flTUl) . f[T2]) . respectively. 

Let us choose Xq, Yq f|T5l) with parameters fixed variationally (see above) as zero approx- 
imation in perturbation theory (fT8|) . (!22|) . By construction of Xq,Yq the emerging pertur- 
bation theory is convergent. Assuming the condition fl26|) fulfilled for the first corrections, 
namely, Ai^x = = v4i, we find the first corrections (pi{C) and gi{r]) as functions of Ai. 
Then we modify the trial function ( |T5|) . 



^So ^ ^o(e)l^J^Hr7) e-^^(«)-^^('') (28) 



and make the variational calculation with this trial function minimizing with respect to 
parameter p. The result is that the optimal value of parameter p remained unchanged with 
respect to the value obtained for the trial function f llSp with extremely high accuracy - within 
10 s.d.! It indicates that the condition fl2B]) is fulfilled with high accuracy. The variational 
energy is changed beyond the 10 s.d. Therefore, the our energies presented in Tables I, II 
are correct in all digits. The separation parameters ,c, Ai-r^ are presented in Table V. It 
allows us to find explicitly Qi{'rf). As an illustration in Figs. 1-4 the functions 

XQ{^)Ylf'\rj) and the first corrections to them are shown for R = 2 a.u. 

Knowledge of wave functions with high local relative accuracy gives us a chance to cal- 
culate matrix elements with controlled relative accuracy < 10^^. As a demonstration we 
calculate the Oscillator Strength as function of interproton distance for the simplest radiative 
transition 2pau — ?■ ^sag are (see e.g. 

/oi(i?) = '^{E'P-^iR) - i?i-«(i?))|Q(i?)|2 , (29) 

where Q(-R) is the matrix element 

Q{R) = (^i^"«(i?)|r|^2pa„^^^^ ^ 

where r is the vector of the electron position measured from the internuclear midpoint, and 
wavefunctions vj/i'^^a^ are given by ( fT5l) . It is assumed this calculation should provide 

at least 5 s.d. correctly. In Table VI the results are presented. For all internuclear distances 
they coincide in 2 s.d. with Bishop et al [7[, thus, indicating the 3rd digit obtained in [7| is 
incorrect for R= 1,2,4 a.u., and in 6 figures with recent results 15^ (with an exception at 
R=l a.u. where it deviates in one unit at the 6th digit) which increases up to 8 figures for 
large R. Modification of (ITSjl by adding the first corrections (|28l) and use it in (1291) does not 
change our 6 s.d. in Table VI. 

Summarizing we want to state that a simple uniform approximation of the eigenfunctions 
for the molecular ion is presented. It allows us to calculate any expectation value or ma- 
trix element with guaranteed accuracy. It manifests the approximate solution of the problem 
of spectra of the molecular ion. In a quite straightforward way similar approximations 
can be constructed for general two-center, one-electron system {Za,Zb,e), in particular, for 
(Hell)^'*". It will be done elsewhere. 

The key element of the procedure is to construct an interpolation between the WKB 
expansion at large distances and perturbation series at small distances for the phase of the 
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wavefunction. Or, in other words, to find an approximate solution for tlie corresponding 
eikonal equation. Separation of variables allowed us to solve this problem. In the case of 
non-separability of variables the WKB expansion of a solution of the eikonal equation can 
not be constructed in unified way, since all depends on the way to approach to infinity. 
However, a reasonable approximation of the first growing terms of the WKB expansion 
seems sufficient to construct the interpolation between large and small distances giving high 
accuracy results. This program was realized for the problem of the hydrogen atom in a 
magnetic field and will be published elsewhere. 

It is worth mentioning a curious fact that the problem ([2]) possesses the hidden algebra 
sl(2) © sl{2). It can be immediately seen - making the gauge rotation of the operators in 
r.h.s. of the equations (E]) and (jH]) with gauge factors e~^^ and e^^, respectively. We obtain 
the operators which are in the universal enveloping algebra of sl{2) (see e.g. 16]). The 
dimension of the representation is —A and —A + ^, respectively. For non-physical values 
of A and integer ratio j the algebras s/(2) appear in the finite-dimensional representation 
realized in action on polynomials in ^, r/. It explains a mystery sometimes observed of the 
existence of polynomial solutions for non-physical values of A in the problem (2) (details 
will be given elsewhere). 
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TABLE I: The total energy Et{R) for Isag state of the -ion compared to [u\ (rounded) 

and Lagrange mesh method. 



R[a.u.] Et[Ry] (Present/ [6] /Mesh) p 



1.0 



L997193 



2.0 



6.0 



10.0 



12.5 



30.0 



40.0 



-0.90357262676 
-0.90357262676 
-0.90357262676 
-1.20526923821 

-1.20526923821 
-1.20526842899 
-1.20526842899 
-1.20526842899 
-1.0239380968 
-1.0239380969 
-1.0239380969 
-1.0011574578 
-1.0011574579 
-1.0011574579 
-1.0002611115 

-1.0002611116 
-1.0000055815 

-1.0000055815 
-1.0000017622 

-1.0000017622 
-1.0000007211 



0.8519936 



1.483403 



1.485015 



3.49506 



5.47987 



6.73221 



15.492 



20.4939 



25.49511 



50.0 



-1.0000007211 
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TABLE II: The total energy Et{R) for 2p(Tu state of the Hg'-ion compared to 6| (rounded) 

and Lagrange mesh method. 



R [a.u.] Et (Present/ [6] /Mesh) [Ry] p 



1.0 



1.997193 



2.0 



4.0 



10.0 



12.54525 



20.0 



30.0 



0.8703727499 
0.8703727498 
0.8703727498 
-0.3332800331 

-0.33328003316 
-0.3350687844 
-0.3350687844 
-0.3350687844 
-0.8911012787 
-0.8911012787 
-0.8911012787 
-0.9998021372 
-0.9998021372 
-0.9998021372 
-1.0001215811 

-1.0001215811 
-1.0000283953 
-1.0000283953 
-1.0000283953 
-1.0000055815 

-1.0000055815 
-1.0000017622 



0.5314196 



1.1536645 



1.155452 



2.3589 



5.47678 



6.75434 



10.4882 



15.492 



20.4939 



40.0 



-1.0000017622 
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TABLE III: The parameters of the function f|T5l) for Isag state of the H^-ion. The 
parameters 0,7 and ai,2)^2,3 are found via minimization. 





i?eo=l-997193 a.u. 


i?=6.0 a.u. 


i?=20.0 a.u. 


a 


1.48407 


3.32381 


10.0453 


P 


1.483403 


3.49506 


10.4882 


7 


1.0299 


0.96357 


0.95774 


ai 


0.9164 


2.597355 


9.8775 


02 


0.05384 


0.53443 


6.8392 


b2 


0.06 


0.588072 


6.9016 


bs 


0.00011 


0.00552 


1.352 



TABLE IV: The parameters of the function (fT5l) for 2pau state of the H^-ion. The 
parameters 0,7 and 01,2,^2,3 are found via minimization. 





R=Q.O a.u. 


Rmin =12.54525 a.u. 


R=20.0 a.u. 


a 


3.24715 


6.5275 


10.7397 


P 


3.43971 


6.75434 


10.4882 


7 


0.95706 


0.97045 


1.03027 


Ol 


2.84566 


6.075 


9.8077 


02 


0.22098 


1.46757 


2.3784 


b2 


0.23611 


1.5349 


2.43705 


b3 


-0.0027 


0.1675 


0.367 
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TABLE V: Separation parameters v4i^^, Ai^^ for lsag,2pau states of the H^-ion compared 

to Scott et al [14]. 





Isag 


2po-„ 


R 


Au 




[14] 


Au 


Ai,n 


[14] 


2.0 


0.811729588 


0.811729585 


0.811729585 


-1.186889395 


-1.186889393 


-1.18688939 


15.0 


48.822353534 


48.822353528 




48.821470973 


48.821470957 




20.0 


90.052891187 


90.052891183 


90.0528912 


90.052877564 


90.052877564 


90.0528776 


30.0 


210.034596601 


210.034596599 




210.034596601 


210.034596599 





TABLE VL Oscillator strength /oi fl29|) for transition 2tot„ — )■ Isag vs R compared to 
Bishop et al p'] and Tsogbayar et al [l5| (rounded). 



R 


present 




[7] 


[15] 




1.0 


0.538675 




0.538 


0.5386739 




1.997193 


0.639595 










2.0 


0.639527 




0.638 


0.6395268 




4.0 


0.469200 




0.476 


0.4692004 




10.0 


2.217 xlO" 


02 


0.022 


2.21706 xlO" 


-02 


15.0 


5.129 xlO- 


04 




5.12939 xlO" 


-04 


20.0 


8.191 xlO" 


06 




8.20513 xlO" 


-06 


30.0 


4.770 xlO" 


09 








40.0 


1.828 xlO" 
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